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Abstract 

QCD running coupling constant a s (m c ) and a s (mb) are determined from 
heavy quarkonia cc and bb decays. The decay rates of V — > 3g and V — ► e + e~ 
for V = J /tp and T are estimated by taking into account both relativistic and 
QCD radiative corrections. The decay amplitudes are derived in the Bethe- 
Salpeter formalism, and the decay rates are estimated by using the meson 
wavefunctions which are obtained with a QCD-inspired inter-quark potential. 
For the V — * 3g decay we find the relativistic correction to be very large and 
to severely suppress the decay rate. Using the experimental values of ratio 
Rg = r^AZZ^rf=j ~ 10) 32 for V = J/ip, T respectively, and the calculated 
widths , we find a s (m c ) = 0.29±0.02 and a s {m^) = 0.20±0.02. These values for 
the QCD running coupling constant are substantially enhanced, as compared 
with the ones obtained without relativistic corrections, and are consistent with 
the QCD scale parameter A^^ w 200MeV. We also find that these results are 
mainly due to kinematic corrections and not sensitive to the dynamical models. 
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1 INTRODUCTION 



Determining the QCD running coupling constant a s at different energy scales is 
very important in the verification of the fundamental theory of strong interaction. 
Among others the heavy quarkonia decays may provide very useful information for 
a s at the heavy quark mass scale. Decay rates of heavy quarkonia in the nonrela- 
tivistic limit with QCD radiative corrections have been studied (see, e.g., refs. [1,2,3]). 
However, the decay rates of many processes are subject to substantial relativictic 
corrections. In particular, the rate of J/ip — > 3g and accordingly the determina- 
tion of a s {m c ) depend rather crucially on the relativistic corrections. Presently the 
ratios of gluonic to leptonic widths of J/ip and T have been precisely measured in 
experiment. However, if using the non-relativistic expressions for decay widths which 
are proportional to |-Rs(0)| 2 , where |-Rs(0)| is the radial wavefunction at the origin, 
and comparing them with the corresponding experimental values, one would get' 2 ' 3 ! 
a s {m L ) = 0.19 and a s (mi>) = 0.17. They are substantially smaller than the expected 
values determined from other experimental results or the QCD scale parameter (for 
a review, see ref.[4]). Theoretically, the difficulty is mainly due to the large rela- 
tivistic effects on the decay widths. It is well known that there are at least two 
important sources of relativistic effects for these processes, one is from the kinemati- 
cal corrections to decay amplitudes, the other is from the bound state wavefunction 
corrections which are concerned with the dynamical effect of quark- ant iquark inter- 
actions. In ref. [5], the first order relativistic corrections were considered based on a 
phenomenological model. But only the relativistic correction originated from kine- 
matics was discussed and no explicit methods with dynamical considerations were 
given to calculate the decay widths. As a result of lack of reliable estimates for rel- 
ativistic corrections to the decay widths in the determination of a s {m c ) and a s (m&), 
either an arbitrary choice for the v 2 /c 2 term is made to get these coupling constants 
enhanced ^1, or the v 2 /c 2 term is arbitrarily neglected but a large effective mass for the 
gluons is introduced' 6 ' . Apparently, a better estimate rather than arbitrary guesses 
for the relativistic corrections is indeed needed, though this is certainly a very difficult 
task without a deep understanding of quark confinement. 

In this paper, as an attempt to tackle this problem, we will use the Bethe-Salpeter 
(BS) formalism to derive the decay amplitudes and to calculate the decay widths of 
V — > e + e _ and V — > 3g, where V is a vector heavy quarkonium state. In this 
approach, the meson is considered as a bound state consist of a pair of constitute 
quark and antiquark, and described by the BS wavefunction which satisfies the BS 
equation. A phenomenological QCD-inspired interquark potential will be used to 
solve for the wavefunctions. These may allow us to give the expressions which take 
into account both relativistic and QCD radiative corrections to next-to-leading order. 
We may then estimate both kinematical and dynamical relativistic corrections to the 
decay widths, and determine the QCD coupling constant from the calculated ratios of 
gluonic to leptonic decay widths. The remainder of this paper is organized as follows. 



2 



In Sec. 2 we study the relativistic corrections to V — > e + e~ and V — > 3g widths. In 
Sec. 3 we discuss the bound state wavefunctions and determine the coupling constants 
a s (m c ) and a s (mj) by using the experimental values of the gluonic and leptonic widths 
of J/ij} and T. A summary and discussion will be given in the last section. 

2 RELATIVISTIC CORRECTIONS TO V — > e+e AND V ^ 3g 

(1). The V — > e + e~ decay 

We first consider the leptonic decay for V — J/tp and T (see also ref.[7]). This 
process proceeds via the QQ annihilation. Define the Bethe-Salpeter wavefunction, 
in general, for a Q\Q 2 bound state \P > with overall mass M and momentum P = 

(Vp 2 + m 2 , P) 

X (xi,x 2 ) =< 0\T^ 1 (x 1 )Mx 2 )\P >, (1) 
where T represents time-order product, and transform it into momentum space 

Xp(q) = e~ tP - x J d 4 xe-^ x x (xux 2 ). (2) 

Here q 1 (mi) is the quark momentum (mass), g 2 (^2) the antiquark momentum (mass), 
and q the relative momentum, 

X = rj 1 x 1 + r] 2 x 2 , x — x 1 — x 2 , (3) 

P = qi + q2, q = mQi - V1Q2, 

where m = ■ r - n ± i — T (i = 1,2). 

In this formalism the quarkonium annilination matrix elements can be written as 

(0 I QIQ I P) = j d'qTr [I(q, P) X p(q)\ , (4) 

where I(q, P) is the interaction vertex of the QQ with other fields (e.g., the photons 
or gluons ) which, in general, may also depend on the variable q°. If I(q, P) is 
independent of q°, Eq. (4) can be written as 

(0 \QIQ\P) = j d 3 qTr[l(q,P)<S> P (q)\, (5) 

where 

M9) = / dq° X p(q) (6) 

is the three dimensional BS wave function of the QQ meson. In the BS formalism 
the wavefunctions in the meson rest frame, where q\ = —q 2 — q, P — (M, 0), can be 
expressed as 

$° P (q) = A^ 7 °(l + 7 )757°A 2 _(-g1^(^) ) 

®p (?) = Ai(g) 7 °(l + l°W^-(-q)f(q), (7) 
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where <3> P ( 1 ) , and $p ( q ) represent the three dimensional wave functions of 0~ and 
mesons respectively, = e^^, e^ is the polarization vector of 1~ meson, tp and / 
are scalar functions which can be obtained by solving the BS equation for _ and 1 _ 
mesons, and A + (A_) are the positive (negative) energy projector operators 

A\(q) = A + (q ± ) = -L(25 + 7 °7 " <?i + ™7°), 
A 2 _(-q) = A_(&) = ^(E-^.q 2 -m 7 °), 



E = y(p + m 2 . (8) 
For process V — > e + e~ with the electron (positron) momentum pi(p2) and helicity 

I(q,P) = -ie^, (9) 
and the amplitude can be written as 

T = e 2 e Q (0 | Q llx Q \ V)u ri { Pl )^v r2 {p 2 )^ (10) 

where M is the meson mass and e Q is the electric charge of the quark Q (Q = c,b). 
Define the decay constant fy by 

j Y Me, = (0 | Q^Q | V) = j d 3 q Tr[^ P (q)}, (11) 

where e^ is the polarization vector of V meson. Then with (7) we can easily find 



2y/3 f 3 ,m + E q 
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^ir/^'— -«pW)' (12) 



where i? = \A? +m 2 , and \/3 is the color factor. Summing over the polarizations of 
the final state and averaging over that of the initial state, it is easy to get the decay 
width 

F(V - e+e-) = ±*a?e 2 Q ft/M. (13) 

Including further the QCD radiative correction and assuming that the radiative cor- 
rection and relativistic correction can be factorized, we then get the following expres- 
sion for the decay width with both relativistic and QCD radiative corrections 

T(V e+O = \*a\\f v lM{\ - l^pl). (14) 
Expanding f v in terms of q /m , to the first order we get 



2. The V -> Sg decay. 

We next study the hadronic decay for V = J/ip and T. We consider a J PC = l~~ 
QQ bound state decaying into three gluons. The decay width is given by 

d<pZ, (16) 



2M 

where the integration / d(f> is over the final-state phase space. Z is defined by 

z = zm 2 , (17) 

where represents summing over the polarizations of the final state and averaging 
over that of the initial state. The decay matrix element T is 

T = -iglTr{T a T b T c ) J d 4 qTrxp(q)I(h, fc 2 , k 3 ; q u q 2 ), (18) 

where 

&(#3 - i 2 + m)fa(4i -fa + ™)h 



I(k 1: k 2: k 3 ;q 1 ,q 2 ) 



((fa ~ Q2) 2 - ™ 2 )((<?i - &i) 2 - m 2 ) 

+ all permutations of 1,2, 3. (19) 



Here fa and e« (i=l,2,3) represent the momenta and polarizations of the three gluons; 
T a ,T b ,T c are the color SU(3) matrices, and a,b,c are the color indices of the three 
gluons; qi and q 2 are the momenta of quark and antiquark, and their time components 
satisfy q® + q® = M. As usual, we assume' 1 ' 5 ' 

q° = qi-q 2 = 0, g° = q° 2 = M/2. (20) 

Thus /(fci, fc 2 , ^3, q±, q 2 ) is independent of q°, and T becomes 

T = -ig 3 s Tr(T a T b T c ) J d 3 qTr{^ P (filfa, k 2 , fa; q)}. (21) 

Substituting (7) and (19) into (21) and (17), we get the expression for Z. In the ex- 
tremely non-relativistic limit the dependence on q of A + (q), A_(g) and I (fa, k 2 , fa; q) 
in Z is neglected and we obtain the usual non-relativistic formula for V — > 3g. In 
order to go beyond the leading order result we expand the matrix element T in terms 
of (f/m 2 . After performing the trace of 7 matrix, for m = M/2 we get to the first 
order of v 2 /c 2 

160<? s 6 f j3 l3 ltA , W1 , ? + q 



J d 3 qd 3 q , {A (x 1 ,x 2 ,x 3 )(l + Q ^ ) 



81M 4 

+A 1 (x 1 ,x 2 , Xz )t±£-}f{^f{^ (22) 
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where f(q) is the scalar wavefunction of vector meson which comes from (7), and 
where 

16 16 32 32 32 16 

A {Xi,X2,Xs) — —J H H ^ H 2~~ 2 

xf ^2 £1X2 xfx2 X1X2 xfa;^ 
+ (two ot/ier permutations) , 

16 16 41 12 12 40 

Ai(xi, X2, X3) — 2 2 ' — 2 ' 2 2~2 

X\ X2 X\X2 X\X2 X1X2 X1X2 

8 8 8 8 48 48 32 

+ ^7 + 



1 2 1 2 x ^ x 9 x ^ x 2 x x 2 x x 2 

+ (two other permutations) , 
x i= ^, i = l,2,3. (23) 

where u;, represents the gluonic energy. The final-state phase space is given by 

After performing the integration, we find the decay width of V — > 3g with relativistic 
corrections to be 

nv 39) = ^p^l c/v , (25) 

where to the first order relativistic correction g v is given by 

= f d 3 q [1- 2.96^] f(q). (26) 

This integral will diverge when the wavefunction / is determined by the inter-quark 
potential which is Coulombic at short distances where one gluon-exchange dominates. 
We could introduce a cut-off scale A of order m{v/c) to regularize the integral. How- 
ever, examining the interaction vertex / given in (19) we can easily see that gy should 
be a convergent quantity if all higher order corrections are taken into consideration, 
because in (19) / becomes inversely proportional to <f as q 1 — > oo. In the absence of 
calculations for higher order corrections we will take the following expression for the 
regulation of gy 

-2 

gv* f d 3 q [1 + 2.96^ I ]- 1 f(q), (27) 
J m A 

which gives the same result as (26) to the first order relativistic correction. Similarly, 
for the leptonic decay the first order expression of decay constant fy in (15) may also 
be regularized in the same manner, which differs from the full expression (12) only 
by higher order terms, and the difference turned out to be small. 
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In connection with (20), another scheme is to take the on-shell condition, which 
assumes (see, e.g., ref.[5]) the quark and antiquark to be on the mass shell 



q\ = q» = M/2 = E = ^m 2 + <f . (28) 

The advantage of this assumption is that the gauge invariance is maintained for the 
on-shell quarks but at the price of treating the quark and antiquark just as free 
particles in a bound state. An apparent problem in this scheme is that if the quark 
mass takes a fixed value (as in the case before) then q 2 will be fixed but not weighted 
by the wavefunction as in the usual bound state description. In order to connect 
the decay process, which occurs at short distances, where quarks are approximately 
on shell, with the bound state wavefunction, which is mainly determined by the 
long distance confinement force, we have to make a compromise between the on-shell 
condition and the bound state description. We will expand the matrix element (21) 
and (17) in terms of — njwv , and allow (f (so the quark mass accordingly) to 
vary in accordance with the bound state wavefunction f(q) which is to be determined 
by the long distance dynamics, or phenomenologically by some dynamical models. 
With this treatment, we get another expression for Z 

160^ 6 r ,o W1 5<f + 5<f 2 , 

+A 1 (x 1 ,x 2 ,x z )^±^)!(q)M), (29) 



where 



16 16 32 32 32 16 

A (xi, x-ii X3) — — H 2 H 2 2 ^ T~2 

X-^ X2 X\%2 X\%2 %1%2 X-^X2 

+ (two other permutations) , 

20 20 28 56 56 88 

Ai(xi, X2, X3) — — 7T H H o o 5—5- 

X\ X2 X\X2 X\X2 X1X2 X1X2 

72 ' 72 16 16 32 



x^X2 x^x<2 x^ X2 x^X2 

+ (two other permutations) . (30) 
Then to the first order relativistic correction, we get 

9v = d 3 q f(q)[l-( Q , - \. +-)- 



2 J 



'8(7r 2 -9) 6 y (M/2) 
^2 

/(?)[! + 2 J^^]" 1 , (31) 



where gy has been regulated in the same manner as (27). 
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Including futher the QCD radiative correction and assuming again that the radia- 
tive and relativistic corrections can be factorized, we then get the following expression 
for the decay width with both relativistic and QCD radiative corrections 

r{V -+ 3g) = 64 ° ( " 2 ^ mo) (1 - C^M, (32) 

where C = 3.7(4.9) for Q = c(6) [2 - 3] . In the extremely non-relativistic limit, the q 2 /m 2 
term is zero, and the scalar wavefunction f(q) satisfies the following relation 



/ d 3 q f(q) 



M 

T #>), (33) 



where ip(0) is the nonrelativistic wavefunction at the origin in coordinate space, we 
then get the well-known zero-order result ' 2 ' 3 ' 5 ' 

nv ^ 3g)= m^p^i mr (34 ) 

With the two expressions (26) and (31), obtained in the two different treatments 
concerning the on-shell condition, we find that they both give very close results. 
This may largely reduce the uncertainties in our calculations associated with the on- 
shell or off-shell descriptions for the quarks which decay at short distances (therefore 
approximately on shell) and are bound together at large distances (therefore off shell). 

Comparing (26) and (31) with (15), we see that the suppression due to relativistic 
correction for V — > 3g is much more severe than for V — > e + e~. This result then rules 
out the conjecture that the relativistic correction to J/ip — > 3g may be neglibibly 
small l 6 '. 



3 DETEMINATION OF a s {m Q ) AND BOUND STATE WAVE FUNC- 
TIONS 

To calculate the widths of leptonic and gluonic decays, we have to know the 
wavefunctions f(q) for cc and bb states, which are determined mainly by the long 
distance inter-quark dynamics. In the absence of a deep understanding for quark 
confinement at present, we will follow a phenomenological approach by using QCD 
inspired inter-quark potentials, which are supported by both lattice QCD calculations 
and heavy quark phenomenology, as the interaction kernel in the BS equation. We 
begin with the bound state BS equation' 8 ' in momentum space 

Hi - m 1 ) X p(qM + m ^ = ^J d ' kG ^ V ~ k )Xp(k), (35) 

where q\ and q2 represent the momenta of quark and antiquark respectively, G(P, q — 
k) is the interaction kernel which dominates the inter-quark dynamics. In solving 
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Eq.(35), in order to avoid the notorious problem due to the excitation of relative time 
variable we have to employ the "instantaneous approximation". Meanwhile, we will 
neglect negative energy projectors in the quark propagators. We then get the reduced 
Salpeter equation' 8 ' for the three dimensional BS wavefunction <&p(q) defined in (6) 

<S> P (q) = p ,_ l Ei _ E / + l J d 3 kG(P, q- £)$p(£) 7 °A 2 „, (36) 

where G(P,q — k) represents the instantaneous potential. We employ the following 
interquark potrntials including a long-ranged confinement potential (Lorentz scalar) 
and a short-ranged one-gluon exchange potential (Lorentz vector) ' 7 ' 

V{r)=V s {r)+ lfl ® 1 »V v (r), 

V s (r) = \r [ - 

ar 

Vv(r) = -l^Me— , (37) 

O T 

where the introduction of the factor e~ ar is to regulate the infrared divergence and 
also to incorporate the color screening effects of the dynamical light quark pairs on 
the QQ linear confinement potential ' 9 '. In momentum space the potentials become' 7 ' 

G(P) = Gs(P)+i»®rG v (P), 
G S (P) = ~-S 3 (P) + } 1 



Gv0) = -±£®-, 08) 
3rr P +a 2 

where a s (p) is the quark-gluon running coupling constant and is assumed to become 
a constant of 0(1) as P ^ 

^ 12 vr 1 
a s (p) = — - 2 . (39) 

11 \n(a + P /A 2 QCD ) 

The constants A, a, a and Aqcd are the parameters that characterize the potential. 
Substituting (7) and (38) into Eq.(36), one derives the equation for the 1~ meson 
wavefunction f(q) in the meson rest frame' 7 ' 

MfM = (E 1 + E 2 )f 1 (q) 

I d 3 k(G s (q- k) - 2G v (q- k))(E im2 + E 2mi )h{k) 

t/lH/O J 



4EiE; 

Ei + E 2 f ,o , „ , _ r*\ E x m 2 + E 2 m x 



J mi + m? 



4E X E 2 J ' mi + m 2 
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= (41) 

The normalization condition of the wavefunction f(q) reads 



where 



/ 



Eq.(40) looks very complicated since higher order terms are all involved. To the 
leading order in the nonrelativistic limit, Eq.(40) is just the ordinary nonrelativistic 
Schrodinger equation. To the first order of v 2 /c 2 , Eq.(40) becomes the well known 
Breit equation with both vector (one-gluon) exchange and scalar (confinement) ex- 
change. 

For the heavy quarkonium cc and bb systems, mi = m 2 = m, Eq.(40) can be 
greatly simplified. By solving Eq.(40) we can find the wave functions for the 1~ 
mesons (see, e.g., ref.[7]). 

Substituting the obtained BS wave functions into (12) and (14), we then get 

r(J/V> — ► e+e~) = 5.QkeV, (43) 

where we have used 

m c = 1.5GeV, A = 0.23GeV 2 , A QCD = 0A8GeV, 

a = O.OQGeV, a = e = 2.7183. (44) 

With these parameters the 2S — IS spacing and J/ip — i] c splitting are required to 
fit the data. Our result is in agreement with the experimental value of T(J/tp — > 
e+e~) = 5.36 ± 0.29keV W. Here in above calculations the value of a s (m c ) in the 
QCD radiative correction factor in (14) is chosen to be 0.29' 3 !, which is also consistent 
with our determination from the ratio of BR(J/ip — > 3g) to BR(J/ip — > e + e _ ) (see 
below). 

Using the experimental data^ 

R > s WJT^H M 10 ' (45) 
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and the calculated widths from (14) and (31), (32), we find 



a s (m c ) = 0.29 ±0.02, (46) 

as campared with the value without relativistic corrections (but with QCD radiative 
corrections) 

a° s (m c ) = 0.19 ±0.02. (47) 

Clearly, it is the strong suppression due to the relativistic correction to J/ip — > 3g 
that enhances the value of a s (m c ) and then makes it consistent with the determined 
QCD scale parameter « 200MeV ( for a review see ref.[4] ). ( Note that a 

slightly larger value for a s {m c ) could be obtained if using (27) rather than (31). ) 

In comparision we have chosen two other groups of parameters and solved the 
BS equation for cc states by requiring again both their 25 — IS spacing and J/ip — 
i] c splitting fitting the data. By the same procedure we get two values of a s (m c ) 
corresponding to the obtained two new wavefunctions of J/ip. With m c = lAGeV, A = 
0.24GeV 2 and other parameters unchanged (the heavy quarkonia mass spectra are not 
sensitive to a and a for a < O.OQGeV) we get a s (m c ) = 0.29. With m c = 1.6GeV, A = 
0.22GeV 2 and other parameters unchanged we get a s {m c ) = 0.28. 

Meanwhile, we have also solved the nonrelativistic Schrodinger equation for the 
scalar function / in the nonrelativistic limit by using the same potentials and parame- 
tres as (38) and (44), and performing the same calculations, we obtain 

a s {m c ) = 0.27. (48) 

Note that in this limit the spin symmetry between the 1~ and 0~ mesons is restored 
and the relativistic correction from the dynamical sourse is eliminated. 

In order to see further the sensitivity of the value of a s (m c ) to the wavefunctions, 
we have also tried the Gaussian function 

f(q) = N exp (49) 

where N is the normalzation factor, and q 2 is the mean value of the momentum 
squared of the quark inside the meson, which may be roughly estimated by using the 
scaling law q 2 . = mC/2 (C = 0.73GeV) found for heavy quarkonia (see, e.g., ref.[10]). 
Then with m c = 1.5GeV we find that while the gluonic and leptonic decay widths 
both become smaller, their ratio does not change very much and gives a s (m c ) = 0.26. 

These results indicate that the determination of a s {m c ) is not sensitive to the 
dynamical corrections to the wavefunctions, and even not very sensitive to the form 
of wavefunctions, and the most important contribution comes from the kinematic 
corrections. 

For the bb system, with a similar calculation for the decay rates of T — > 3g and 
T — > e + e _ with both relativistic and QCD radiative corrections taken into account, 
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and using the observed value of the ratio' 4 ' 



we find 

a s {m b ) = 0.20 ±0.02, (51) 

as compared with the value without relativistic corrections (but with QCD radiative 
corrections) 

a° s (m b ) = 0.17 ±0.02. (52) 

The value given in (51) agrees with the expected value from the QCD scale parameter. 
Moreover, with both relativistic and QCD radiative corrections, and a s (m b ) = 0.20, 
we get' 7 ' 

r(T -> e + e~) = (1.3±0.2)tfeV, (53) 

which is also in agreement with data' 4 '. 

We may then conclude that by estimating the relativistic corrections to the gluonic 
and leptonic decays of heavy quarkonia, we find that the relativistic effects substan- 
tially suppress the V — > 3g decays, and consequently the determined values of the 
QCD running coupling constant at the heavy quark mass scale can get enhanced, and 
can be consistent with other theoretical and experimental expectations. 



4 SUMMARY AND DISCUSSION 



In this paper we calculated the first order relativistic corrections to V — > e + e _ and 
V — > 3g for the vector heavy quarkonia V = J/ip and T, based on the BS formalism 
for the decay amplitudes and bound state wavefunctions. We derived the coefficients 
of the v 2 /c 2 term in the decay rates with two different treatments for the on-shell con- 
dition of bound quarks, and obtained very similar values for the coefficients. These 
results may largely reduce the uncertaity in our calculations concerning the on-shell 
condition in these decay processes. To maintain gauge invariance in the decay ampli- 
tudes, the on-shell quarks are certainly better than off-shell quarks. But within the 
quark- ant iquark sector of the Fock states for a heavy quarkonium, gauge invariance 
may not be guaranteed and other higher Fock states may be needed. Nevertheless, 
the closeness of the two large negatrive values of the coefficients of the (f/m 2 term 
(see (26) and (31)) in the V — > 3g widths within the two treatments may encour- 
age us to conclude that the relativistic effects suppress V — > 3g decays much more 
severely than V — > e + e _ decays, therefore can make the coupling constants a s {m c ) 
and a s {m b ) substantially enhanced, as compared with the values obtained without 
relativistic corrections. 

For a more accurate estimate of relativistic corrections, higher orders, e.g., the 
(v 2 /c 2 ) 2 terms should be taken into account, but this is difficult in our approach, 
since to the higher orders many other effects e.g. the quark pairs and the dynamical 
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gluons may be involved. In connection with higher order effects, we have regularized 
the singularity associated with the derivative of the wavefunctions at origin in a 
simple manner, which is valid to the first order of v 2 /c 2 but is uncertain to higher 
order corrections. 

We have solved the BS equation for the bound state wavefunctions with QCD 
inspired inter-quark potentials (linear confinement potential plus one gluon exchange 
potential) as the BS kernel. With some popular parameters for the potentials we ob- 
tained the wavefunctions and used them to calculate the gluonic and leptonic decay 
widths and their ratios. By comparing the BS wavefunctions with Schrodinger wave- 
functions and Gaussian-type wavefunctions, we found that the ratios are not sensitive 
to the dynamical relativistic effects on the wavefunctions, and even not very sensitive 
to the form of wavefunctions. We may then conclude that the ratios are insensitive to 
the dynamical models, and the relativistic effects on the ratios are mainly originated 
from the kinematic parts of decay amplitudes. 

As for the numerical result, using the experimental values of ratio R g = fT^ zz^ F ) 
10, 32 for V = J/ip, T respectively, and the calculated widths, we found a s (m c ) = 
0.29 ± 0.02 and a s {m b ) = 0.20 ± 0.02. These values for the QCD running coupling 
constant are consistent with the QCD scale parameter pa 200MeV. 

Recently, there have been some significant progresses in the study of heavy quarko- 
nium decays based on a more fundamental approach of NRQCD (nonrelativistic QCD) 
(see refs. [11,12]). Many important issues were clarified by this study. It will be inter- 
esting to compare their results with ours in connection with the gluonic and leptonic 
decay widths and the determination of the strong coupling constant at the heavy 
quark mass scales from these decays. 
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